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ABSTRACT 


These  investigations  generalize  Miyata's  synthesis  of  passive 
driving-point  impedance  functions  by  developing  a  step-by-step  com- 
putational technique  for  augmenting  a  general  driving-point  impedance 
to  insure  that  the  real  part  of  the  augmented  impedance  is  positive 
term  by  term.   This  goal  has  been  accomplished  and  programmed  under 
the  provision  that  the  real  part  of  the  original  impedance  is  minimum 
reactive,  i.e.,  has  no  zeros  on  the  joo-axis.   The  latter  requirement 
is  not  restrictive  since  zeros  on  the  joo-axis  can  be  removed  a  priori. 
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I.   INTRODUCTION 

7  8 
The  Miyata  '   method  of  transformerless  synthesis  provides  a  method 

for  separating  a  given  impedance  into  additive  components  which  indi- 
vidually satisfy  certain  special  conditions  that  permit  their  realization 
in  simple  terms.   This  method  of  separation  is  based  upon  the  even  part 
of  the  given  impedance  and  can  be  stated  briefly  as  follows :   Given  a 
minimum  reactive  driving  point  impedance  Z(s),  find  R(oo)  the  real  part 
of  Z(s).   Then  split  R(cd)  into  a  sum  of  even  functions  R  (cd)  ,  (p=0 , 1 ,  .  .  .  ,n) 
each  satisfying  the  requirements  of  ReZ  (jcd)  .   Find  Z  (s)  corresponding  to 

each  R  (cd)  .   The  realization  of  Z(s)  reduces  to  the  realization  of  each  of 
p 

the  Z  (s)  and  the  series  connection  of  the  resulting  one  ports  completes 
the  synthesis. 

An  inherent  difficulty  in  the  above  technique  lies  in  the  require- 
ment that  the  R  (to)  be  positive  for  all  cd.   To  overcome  this  difficulty 
Miyata  suggested  the  use  of  surplus  factors.   An  example  will  illustrate 
this  method. 

Let  Z(s)  be  given  as 

z(s)=6sW9  (11) 

s  +3s+4 

from  which  we  obtain 

4   2 
R(o3)  =  ^  -to2+36  =  R2(cd)  +  r  (a>)  +  RQ(o))         (1.2) 
B(co  ) 

This  so  called  "single-n  split"  method  cannot  be  applied  because 

R,  (cd)  is  negative.   Therefore  use  a  surplus  factor  (s+1)  to  obtain  an 

augmented  impedance,  Z  (s). 

7  /.\  ■   s+1  vr    ^  -  6s3+15s2+18s+9  ,.    „. 

Za(s)  "  i+1  Z<s>  =  3^  2M1  Z  (1,3) 

s  +4  s  +7s+4 


for  which  the  numerator  of  ReZ  (jo>)  contains  only  positive  terms. 

8 

-or   \        6cd  +30ao  +36  . 

R(od)  =  j (1.4) 

B'(od  ) 
Obviously  the  use  of  surplus  factors  complicates  matters  since 

higher-order  factors  are  often  required. 

1  8 
At  this  point  most  authors    considering  the  technique,  including 

Miyata,  generally  abandon  the  use  of  surplus  factors  without  discussing 
any  general  method  for  their  selection  and  proceed  to  describe  various 
groupings  of  the  individual  positive  and  negative  R  (cd)  which  remain 
positive  for  all  <x>.      These  grouping  methods  generally  require  factori- 
zation of  the  numerator  with  its  attendant  difficulties  and  no  general 
method  is  available  in  the  literature  nor  is  it  clear  that  such  groupings 
are  always  possible. 

The  purpose  of  this  thesis  is  to  investigate  in  detail  the  use  of 
surplus  factors  to  develop  augmented  impedances  that  are  positive  term 
by  term. 

Chapter  two  is  a  general  consideration  of  the  Miyata  method  which 

describes  the  reasoning  behind  the  technique  and  the  mechanics  of 

developing  the  electric  circuit  and  component  values  from  the  original 

driving-point  impedance  when  all  terms  in  the  numerator  of  ReZ(ja))  are 

initially  positive.   The  significant  results  of  this  section  include: 

The  circuit  realization  is  transformerless,  an  important  advantage 
over  the  Bott-Duffin  technique. 

The  complete  circuit  development  consists  of  a  simple  unbalanced 
ladder  network  with  series  L's  and  shunt  C's  terminated  in  a  resistance. 

By  term  grouping  without  factorization,  the  number  of  circuit  elements 

required  is  less  than  for  the  modified  Bott-Duffin  for  all  driving-point 

(dp)  impedance  functions  of  second  order  or  greater.   The  saving  in 
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number  of  elements  becomes  quite  dramatic  for  higher-order  polynomials 
so  that  even  though  the  order  of  the  original  dp  impedance  must  often 
be  increased  by  the  use  of  surplus  factors,  large  savings  are  available. 

In  general,  chapter  two  is  included  to  provide  sufficient  incentive 
to  undertake  the  task  of  developing  a  general  method  of  selecting 
surplus  factors  of  minimum  order. 

In  chapter  three,  starting  with  a  dp  impedance  that  is  the  ratio  of 
two  general  polynomials  of  the  same  order  in  s  and  a  general  augmenting 
polynomial  (surplus  factor),  the  real  part  of  Z(s)  is  developed.   Since 
the  resultant  coefficients  of  od   must  be  non-negative,  these  coefficients 
are  the  required  constraints  that  must  be  satisfied  to  insure  that  ReZ(jco) 
is  negative  term  by  term.   These  coefficient  constraints  constitute  a 
set  of  non-linear  inequalities  in  the  variables  6.  (i=0, 1, . . . ,n)  which 
are  the  coefficients  of  the  augmenting  polynomial  (i.e.,  P(s)=6  +6.S+...+6  s  ) 

Also  in  this  section,  the  non-linear  inequalities  developed  are  linearized 

9  2 

using  a  set  of  congruent   transformations  and  solved  by  the  simplex  method. 

The  difficulties  involved  in  the  inverse  transformation  are  circumvented  by 
considering  the  form  of  the  solved  constraints  which  are  the  coefficients 
of  the  real  part  of  the  augmented  dp  impedance.   The  process  is  completed 
by  demonstrating  how  a  knowledge  of  the  form  of  the  solved  constraints 
allows  one  to  find  the  numerator  of  the  real  part  by  simple  long  division. 

■ 

Determination  of  the  denominator  is  straightforward  and  the  problem  is 
solved . 

A  digital  computer  program  for  part  of  the  routine  described  above  is 
developed  in  chapter  four,  and  a  sample  problem  worked  in  section  five. 

In  summary,  a  method  has  been  developed  which  will  invariable  find 
an  optimal  augmenting  polynomial  such  that  the  resultant  augmented  dp 


impedance  has  a  real  part  which  is  positive  term  by  term.   The  method 
is  optimal  in  the  sense  that  as  many  coefficients  as  possible  of  the 
real  part  are  zero,  thereby  reducing  the  required  number  of  circuit 
elements  to  a  minimum. 

The  value  of  this  method  lies  in  the  fact  that  it  will  invariably 
produce  a  positive  real  part  for  the  augmented  impedance  in  a  step-by- 
step  fashion.   In  fact,  the  entire  process  could  be  programmed  such  that 
for  any  input  minimum-reactive  dp  impedance  the  output  would  be  the 
required  circuit  with  accompanying  element  values. 

The  single  major  disadvantage  is  that  the  augmented  dp  impedance  may 
have  a  large  number  of  terms  which  would  consequently  require  a  large 
number  of  circuit  elements  to  realize  (the  circumstances  that  determine 
this  possibility  are  discussed  in  chapter  two). 
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II.   GENERAL  CONSIDERATION  OF  THE  MIYATA  SYNTHESIS  TECHNIQUE 


The  Miyata  method  for  realization  of  a  passive,  minimum-reactive, 
driving-point  impedance  proceeds  as  follows :  '   For  an  impedance  having 
no  joD-axis  zeros  or  poles  written 

m1  +  n 

Z(s)  =    .  n  (2.1) 

m2  +  n2 

where  m  denotes  an  even  polynomial  in  s  and  n  denotes  an  odd  polynomial 


in  s 


Its   real   part    for  s=jco  is 


m..m_    -   n,n_ 

ReZ(jco)   =-±y Y^  (2-2) 

nu      -   n_ 


or 


aia2                   2n  .             .2  A2n 

A  +A..CO  +.  .  .+A  co  A            A  od  A  en 

nn/.\           ol                    n  o.i           .  .n                    /oo\ 

ReZ(ja)) J 27T  =  T  +   2~  +  *  '  +  T  (2-3) 

B  +B.0D  +.  .  .+B  cd  B(co  )    B(co  )  B  (cd  ) 
o   l        n 

Each  of  the  terms  on  the  right-hand  side,  when  the  coefficients 
A  ,A. , . . . ,A  are  positive,  may  be  shown  to  generate  an  impedance  that 
can  be  realized  without  mutual  inductive  coupling.   The  series  con- 
nections of  such  separate  realizations  yields  the  desired  synthesis  of 
Z(s). 

To  see  why  this  is  so,  consider  the  component  impedances  of  the  real 
part.   The  first  term  has  n  double-order  zeros  at  s=  infinity;  the  second 
has  one  double-order  zero  at  s=0  and  n-1  zeros  at  s=  infinity;  the  third 
term  has  two  double-order  zeros  at  s=0  and  n-2  zeros  at  s=  infinity; 
finally,  the  last  term  has  all  its  n  double-order  zeros  at  s=0.   There- 
fore, all  of  these  terms  have  their  double-order  zeros  located  exclusively 
at  s=0  and  s=  infinity. 
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When  the  real  part  of  an  impedance  is  zero  at  some  jco-axis  point, 
it  need  not  necessarily  follow  that  the  entire  impedance  is  zero. 
However,  the  joo-axis  points  under  consideration  are  restricted  to  s=0 
and  s=  infinity  and,  since  a  minimum-reactive  function  must  be  purely 
real  at  these  points,  it  follows  that  when  the  real  part  is  zero  the 
whole  impedance  must  necessarily  be  zero. 

Next  we  use  the  Darlington   theory  to  show  that  these  impedances  are 
completely  developable  in  simple  terms.   The  Darlington  theory  shows  that 
in  a  cascade  development  of  an  impedance  each  component  lossless  network 
produces  one  set  of  zeros  out  of  the  total  complement  of  zeros  which  the 
even  part  of  that  impedance  contains.   Therefore,  the  network  generated 
by  the  first  cycle  provides  one  set  of  zeros  and  leaves  a  remainder  whose 
even  part  contains  the  remaining  even-part  zeros.   The  next  network  in 
the  chain  provides  another  set  and  leaves  a  second  remainder  with  the 
leftover  even-part  zeros  -  and  so  forth. 

Applying  this  reasoning  to  the  real  part  of  the  first  term  of  the 
expansion  of  the  real  past  of  Z(s),  we  can  predict  that  the  corres- 
ponding impedances  will  have  a  zero  at  s=  infinity  and  its  reciprocal 
a  pole  at  this  point.   Removal  of  this  pole  as  a  shunt  capacitor  is  one 
Darlington  cycle.   The  admittance  remainder  after  this  pole  at  s=  infinity 
is  removed  has  an  even  part  that  still  has  n-1  double-order  zeros  at 
s=  infinity  since  so  far  only  one  has  been  used  up.   This  admittance 
remainder  must  have  a  zero  at  s=  infinity  and  its  reciprocal  must  have 
a  pole  there.   Removal  of  this  pole  as  a  series  inductor  constitutes  the 
second  cycle.   The  remainder  now  has  an  even  part  with  n-2  double-order 
zeros  at  s=  infinity  and  is  zero  at  s=  infinity.   Its  reciprocal  must 
have  a  pole  at  this  point  and  its  removal  as  a  shunt  capacitance  completes 
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the  third  cycle.   This  process  continues  yielding  series  inductive  and 
shunt  capacitive  branches  of  total  number  n,  the  number  of  double-order 
zeros  possessed  by  the  original  real  part.   When  the  last  of  these  zeros 
has  been  accounted  for,  the  remainder  is  reduced  to  a  constant  because 
its  real  part  has  no  more  zeros.   The  complete  development  consists  of  a 
simple  unbalanced  ladder  network  with  series  inductances  and  shunt  cap- 
acitances terminated  in  a  resistance.   Subsequent  terms  are  similarly 
developed. 

Using  this  single-n  split  technique,  the  total  number  of  reactive 

2 

elements  is  n  ,  where  n  equals  the  highest  degree  of  either  the  numerator 

or  denominator  polynomial  in  s.   Comparison  with  the  number  of  elements 

required  for  a  Bott-Duffin  synthesis  is  displayed  in  Table  I  along  with 

the  reductions  due  to  the  1/2-n  split  technique  which  is  subsequently 

discussed.   The  geometric  character  of  the  Bott-Duffin  method  compared 

with  the  arithmetic  character  of  the  Miyata  method  is  strongly  apparent 

in  this  table.   Even  though  special  devices  are  required  to  salvage  the 

Miyata  method  when  all  the  terms  of  the  real  part  numerator  are  not 

initially  positive,  there  appears  to  be  considerable  advantage  to 

exploring  the  Miyata  procedure  where  more  elaborate  impedance  functions 

are  involved.   Before  one  gets  too  enthusiastic  about  the  Miyata  method, 

it  would  be  well  to  consider  the  requirement  that  all  the  coefficients 

in  the  polynomial  formed  by  the  numerator  of  the  real  part  of  the 

driving-point  impedance  must  be  positive. 

A  sufficient  condition  that  the  coefficients  be  positive  can  be 

obtained  by  recalling  that  a  Hurwitz  polynomial  has  positive  coefficients 

Thus  if 

A(od2)  =  An  +  ...  +  A  u)2n  (2.4) 

0  n 
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2 

has  no  zeros  in  the  right-half  od  plane,  then  A  ,  .  .  .  ,A  are  positive 

0  n 


Or   if 


2 
A(-s    )   =  nunu    -   n1n2  (2.5) 


2 

has   no  zeros    in   the    left-half   s      plane,    then   the  A    ,  .  .  .  ,A     are    positive. 

2 

Since    the    left-half   s      plane  maps    into  a  region   of   the   s    plane    con- 
taining all   points   not  more   than  plus    or  minus   45      from  the    jco-axis  ,   we 
can   say  the   coefficients    in  eq.    2.4  are    positive    so    long  as    there   are   no 
even-part   zeros    of  Z(s)    closer   than  45      to   the   jco-axis  . 

Although   the   coefficients    in  2.4  may   still  be    positive  when   some 
even-part   zeros    lie   closer   to   the    joo-axis    than   this    sufficiency  con- 
dition  permits,    the   chances    of   this   desired   condition  being   fulfilled 
when  even-part   zeros   are   relatively  close   to  the   jco-axis   become   remote 
and  vanish   completely  when  a   single    pair   of  even-part   zeros    lie   upon   the 
ju)-axis  . 

In  a  manner   of  speaking,   we   can   say  that   a   polynomial   that   has    some 
zeros   above    the   45      line    in   the   s    plane  may  still  yield  all   positive 
coefficients    if   there   are   a   sufficient   number   of  zeros   below  the   45 
lines    to   prodice   a   compensating  effect.      Therefore,    in  a   situation    in 
which   the   even   part    of  Z(s)   has    some   negative   coefficients,   we   can  augment 
the   real   part  with   appropriately   chosen   factors    to  bring  about   the   desired 
compensation. 

The  method   of   compensation   fails    if   the   even   part    of  Z(s)   has   any 
joo-axis   zeros,    since   the   number   of  necessary  compensating   factors    tends 
toward    infinity  as   one   or  more   pairs    of   even-part    zeros   approach   the 
joo-axis.      With   such  a   compensating   scheme    the   number   of   terms    in   the   real 
part   and  hence    the   number   of  necessary  elements    increases   rapidly. 
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One   method   to  reduce    the   number   of  elements    is   considered   now,    the 

1   8 
1/2-n  split.    '        Miyata   points    out    that   the   detailed   procedure    for 

carrying   out   the   basic  method   of  decomposition  may  be   varied   by   separating 

the   real-part   numerator   polynomial    into  groups    of   terms   rather   than 

separate    terms. 

Consider 

/a    i  a      2,         ,.           2(n-k)N         ,.               2 (n-k+1),                  2nN 
(A  +A_oo  +.  .  .+A     ,oo            ')         (A            oo                   +  .  .  .+A  oo      ) 
,-.   r, , .    x             o      1                   n-k                      .         n-k+1                                   n  ,n    r  \ 

ReZ(joo)   =  5 +  s (2.6) 

B(oo  )  B((D  ) 

where  k=n/2  (n  even)   k=(n+l)/2  (n  odd). 
S.  N.  Hunt  ,  who  refers  to  this  decomposition  as  the  1/2-n  split,  finds 
that  it  yields  a  smaller  total  number  of  elements  than  any  other  division 
of  the  numerator  polynomial  into  parts.   It  requires  fewer  elements  than 
the  Bott-Duffin  method  for  any  degree  n  equal  to  or  greater  than  2.   The 
impedance  generated  by  the  first  term  is  developable  into  a  simple  unbalance 
ladder  with  series  L's  and  shunt  C's  (total  number  k)  terminated  in  a 
remainder  whose  even-part  zeros  are  the  same  finite  non-zero  ones  contained 
in  the  even  part  of  this  impedance.   The  second  term  is  developable  into  an 
unbalanced  ladder  with  series  C's  and  shunt  L's  (total  number  n-k+1)  and  a 
terminal  impedance  whose  even-part  zeros  are  those  that  are  contained  in 
the  even  part  of  the  second  impedance.   To  the  terminal  impedances  which 
appear  as  remainder  functions  we  can  apply  a  similar  treatment  and  continue 
until  the  ultimate  remainder  functions  are  resistances. 

The  process  of  augmentation  by  surplus  factors  is  also  necessary  in 

order  that  Miyata 's  elegant  method  of  determining  an  impedance  from  its 

7  8 
real  part  may  be  used.   The  method  proceeds  as  follows:  ' 

m1  +  n1 

Z(s)  ■  .   ■  -1  (2-7) 

m„  t  n_ 
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and 


mm   -  n  n    M/  * 

M(.)  -  -V— p  -  fgi  .  (2.8) 

m2  "  n2 


Next  consider  an  auxiliary  function 


N(m  '  +  n  •) 

Z'(s)  =  .  (2.9) 

m_  t  n„ 


Its   even  part    is 


N(m, 'nu    -   n    'n~) 

EvZ'(s)   ' ~-= 2  .  (2.10 

iru      -    n_ 

Set  the  two  even  parts  to  be  equal; 

EvZ'(s)  =  EvZ(s)  (2.11) 

which  can  be  done  provided 

m  'nu  -  n.'n  =  1.  (2.12) 

where  the  degree  of  m, ,+n1 '  is  that  of  m1+n1 .   Polynomials  m, '  and 
n, '  can  always  be  found  to  satisfy  this  requirements.   Z'(s)  however  will 
have  a  numerator  of  degree  two  or  more  higher  than  its  denominator  and  is 
therefore  not  positive  real  (p.r.).   To  find  the  p.r.  function  Z(s),  carry 
out  a  long  division  on  Z'(s)  until  the  degree  of  the  numerator  does  not 
exceed  the  degree  of  the  denominator.   That  is, 

N(m  '  +  n  ')  m  +  n 

Z'<S)  =   m  +  n      =  q(5)  +  J  +   nr  (2'13) 

m2  +  n2  m2  +  n2 

where  q(s)  is  the  quotient  polynomial  and  m  +n   is  the  remainder. 


Let  Z  (s)  be 
r 


so  that 


m  +  n 
r      m„  +  n„ 


Z(s)  =  Zr(s)  +  q(s)  (2.15) 
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EvZ(s)  =  EvZ  (s)  +  Evq(s) 


Recalling  that  EvZ'(s)  =  EvZ(s), 

(2.16) 

If  the  polynomial  q(s)  contains  an  even  power  of  s  then  the  right 
hand  side  of  eq.  2.16  tends  to  infinity  as  s=jao  tends  to  infinity.   On 
the  other  hand,  EvZ(s)   .   =  R(cu)  is  bounded  at  infinity.   Therefore  q(s) 

S=J0D 


must  be  an  odd  polynomial  and  as  a  result 


EvZ  (s)  =  EvZ(s) 


(2.17) 


Since   the   numerator   of  Z    (s)   does   not   exceed   the   denominator   Z    (s) 
has   no   pole  at    infinity   or  at   any   other   finite    frequency   on   the   jco-axis 
since    its   denominator   is    that    of  Z(s). 

Briefly  Z    (s)   has   no   poles    in   the   right-half   plane   or   on   the    jto-axis 
and    its   even   part    for   s=joo  is   always   non-negative.      Therefore   Z    (s)    is 
positive   real.      To   illustrate: 


Given  ReZ    (joo)   = 


1  +  OD 


(2.17a) 


7    (   \   -  2/3s     +  4/3s  +   1 
Z    (s)   -  -= j 

s     +  2s     +  2s  +   1 


(2.17b) 


from  any  method  whatsoever.      Suppose  we  want   an   impedance  with   the 

real   part   ReZ-Qoo)    =  go  /1-kD     which    is   2.17a  multiplied  by  a)     and    implies 

2 

m=-s    . 

2 

Then  F(s)   equals   eq.    2.17b  multiplied   by   -s      and   the    long   division 

takes    the    form 


3  2 

s     +  2s     +  s  +   1 


and 


-2/3: 


■2/3s,    -   4/3s_    -      s      9 
■2/3s      -   4/3s      -   4/3s 


l/3s     +  2/3s 


l/3s^  +  2/3s 
L\KS)  3  2 

s     +  2s     +  2s  +   1 
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Using  this  procedure  any  Z  (s)  =00  /(1+od  )  n=l,2,...  can  be 

found  from  Z  (s).   If  we  want  an  impedance  with  the  real  part 

2      4      6 
ReZ(jo))  =  A  +  A-cd  +  A.cu  +  A0cd 
VJ  '    o    1      2      3 

Z(s)  may  be  generated  as 

A  Z  (s)  +  A  Z.Cs)  +  A  Z0(s)  +  A_Z.(s) 
00       ii       ll.  Jj 


where  continuing  as  above 


7  ,  .        2/3s2  +  l/3s 
Z0(s)  = 


2X  *  3     2 

s  +  2s  +  2s  +1 


„  ,  N   s3  +  4/3s2  +  2/3s 
Z^(s)  =  —z x 

s  +  2s  +  2s  +  1 


2 
Z_(s)  is  found  by  multiplying  Z,(s)  by  -s   and  Z~(s)  found  by 

2 

multiplying  Z„(s)  by  -s   and  performing  the  indicated  long  division 
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III.   GENERAL  STEP-BY-STEP  PROCEDURE  FOR  THE  SELECTION  OF  SURPLUS  FACTORS 


The  topic  of  this  chapter  is  the  step-by-step  procedure  for  develop- 
ing the  real  part  of  an  augmented  impedance  that  is  positive  term  by  term 
from  a  general  n   -order,  passive,  minimum-reactive,  driving-point 
impedance  augmented  by  a  general  m   -order  augmenting  polynomial. 


A.   CONSTRAINT  DEVELOPMENT 

This  section  is  devote'1  to  the  development  of  the  required  constraints 
to  insure  a  positive  real  part  for  the  dp  impedance  and  begins  with  a 
general  impedance  function  and  augmenting  polynomial  as  follows: 

Consider  first  the  general  driving  point  impedance  Z(s) 

N 

E      q 

a  s^ 


{    }  N  Q   q   D(s)  ^*i; 


q=0 
The  general  augmenting  polynomial  P(s) 

T 


P(s)  =  E  6  sr  (3.2) 

r=0 


The  product  of  Z(s)  and  P(s)/P(s)  is  the  augmented  driving  point 


impedance  Z  (s)  ; 

3 


N        T 

S  (a  sq  E  6  sr) 

i:  (3  sq    e    6rsr) 

q=0  q   r=0 
N 
E  a  6rsq+r 

T 

E  3  6  -q+r 


q=0  "I'*8 
r=0 
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The  even  part  of  the  augmented  driving  point  impedance  is 


mlm2   -"l"2   .  N^s2 
EvZa(s)  "   ,2     ,2     -D'(s) 

2    "   2 


(3.5) 


where 


(N+T)/2 

T  9 

mn   =  S  CU    6  s 
1     __   2u-r  r 
u=0 

r=0 


(N+T-l)/2 
T 

1      _   2u+l-r  r 
u=0 

r=0 


2u+l 


(N+T)/2 

T 

m,'  =  E  g0   6s 
2      _   2u-r  r 
u=0 

r=0 


2u 


(NfT-l)/2 

T 

n  '  =  S  0      6  s 
2     rt  K2u+l-r  r 
u=0 

r=0 


(3.6) 
2u+l 


and  noninteger  values  of  (IW-T)/2  and  (N+T-l)/2  are  rounded  to  the 
next  lowest  integer. 


Expanding  the  even  part  of  Z  (s) 

3. 


(N+T)/2 
T 


(N+T)/2 
T 


(N+T-l)/2 
T 


(NT-l)/2 
T 


EvZ   (s)  =     S  cu       6   s2u     S  3_       6   s2u  -     2  a,    , .      6   s2u+1  Eg,  6s 

ax  .   2u-r  r  _^2u-r  r  _  2u+l-r  r  _^2u+l-r  r 

u=0  u=0  u=0  u=0 


2u+l 


r=0 


r=0 


r=0 


r=0 


D'(s) 


(3.7) 


Combining  products, 

(W-T)/2  (W-T-l)/2  (3.8) 

EvZ   (s)  =      I     s2(u+v)a9       39       6  6     -       I     s2(u+v+1)a9^     39  +1     6  6 
3  u,v=0  2v_r   2u"q  r   q        u,v=0  2vfl-rp2u+l-q   r   q 

q,r=0  q,r=0 

D'(s) 


Combining  the  difference 


(N+T)/2 
T 

EvZ   (s)  =      E     s2(u+v)6  6(ct9       09       -s2a9    ,  .     0,..      ) 
a  rq2v-r2u-q  2v+l-r   2u+l-q 

u,v=0  n  — _7.    N M 

a  D'(s) 

r,q=0 


(3.9) 
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Since  ReZ  (jo))  =  EvZ  (s)  _.   and  the  Miyata  method  requires  that 
a  *        a   s—jco 

the  real  part  be  positive  term  by  term,  constraints  must  be  developed 
to  insure  this  requirement.   First,  consider  the  denominator  D'(s)  above 

D'(s)  =  m2'2  -  n2'2  =  D(jo))D(-jo))  =  |d(jcd)|  2  (3.10) 

which  is  never  negative  and  is  therefore  of  no  concern. 

Thus  the  constraint  is  reduced  to 

N'(tD2)    =  m    'm    *    -   n/n    '  =0   for   O^iio  (3.11) 

12  12      s=ja) 

The   coefficient    of  each   term   in  eq.    3.9  must   be   non-negative   which 

requires  that 

(N+T-l)/2 
(NfT)/2 

E  s    6  6  (a,  .„  0   39    -  a,  .,  ,   P0  .,  )  =  0     (3.12) 

r  q  4m+2-2u-r  2u-q    4m+l-2u-r  2u+l-q 
m=0 

u=0 

r,q=0 

and 

(N+T-2)/2 

(W-T)/2 

T 

4m+4  > 

2  s    6  6  (a,  .,  0      30    -  a,  .«  0   30  ,,  )  =  0    (3.13) 

_       r  q  4m+4-2u-r  2u-q    4m+3-2u-rK2u+l-q 
n=0         -ii  -i 

u=0 
r,q=0 

In  equation  3.9  when  the  exponent  2(u+v)  =  4m+2  (m=0, 1, . . . , (NfT- l)/2) , 
the  summation  constraint  becomes  eq.  3.12;  and  when  the  exponent  2(u+v)  = 
4m+4  (m=0, 1,  .  .  .  ,  (IW-T-2)/2)  ,  the  summation  constraint  becomes  eq.  3.13 
since  the  summation  must  be  non-negative  when  evaluated  at  s=jco. 

Collecting  terms  and  simplifying,  eqs.  3.12  and  3.13  become 

)  (a  ,B  , ,  -  a 
q  y-3  €+1    7 


E     6  Ma  i  ,,  -  a  i  )  =  0  (3.14) 

u,r,q,m,-0 


E    6  6  (a  3  -  a  ,p  ,.)  =  0  (3.15) 

_   r  qN  ywP,  7-I  e+17 

u,r,q,m=0    n  '  '        ° 
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where 

u=0,l, . . . ,(N+T-l)/2  u=0,l, . . . ,NfT/2 

r=0,l,...,T  r=0,l,...,T 

q=0,l,...,T  q=0,l,...,T 

m=0,l,. . . ,N+T-l/2  m=0,l,. . . , (ttfT-2)/2 

(3.14)  (3.15) 

and 

7  =  4m+4-2u-4    £   =  2u-q 
N  =  Order  of  original  impedance.   T  =  order  of  augmenting  polynomial. 
The  final  value  of  the  indices  is  chosen  so  that  if  the  value  is  non- 
integer  the  rule  is  to  round  down  to  the  next  lower  integer. 

In  the  summations  above,  given  values  for  a  and  |3 ,  6's  must  be  chosen 
so  that  these  inequalities  are  satisfied. 

Expanding  the  generalized  constraints  for  a  first-order  augmenting 
polynomial  (i.e.,  T  =  1;  r,q  =  0,1) 


where 


A6Q2  +  £6^  +  C6X2   =  0 
D6Q     +  Eo^  +  lbx2  =  0 


-  ==  a4m+l-2uP2u+l   "   a4m+2-2uP2u 


-  M  a4m-2ue2u+l   "   a4m+2-2u32u-l 


-  I=  a4m-2ue2u  "a4m+l-2u32u-l 


-  ==  a4m+4-2ue2u  "  a4m+3-2u32u+l 


£  zz   a4m+4-2ue2u-l  '  a4m+2-2u32u+l 


-  ==  a4m+3-2u32u-l  "  a4m+2-2u32u 


(3.16) 


(3.17) 
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A  = 


alPl   "   CX230   "  a082 


a333  +  aL35  +  a331  -  a234  -  a4P2  -  aQ36  -  a63Q 


B  =  0 


(3.18) 
(3.19) 


C   = 


a4P0  +  a034  +  a232    "  a33l   ■   aiP3 


(3.20) 


D  = 


E  =  0 


(3.20a) 


(3.20b) 


F  = 


(3.20c) 


and  so  for  a  first-order  augmenting  polynomial  the  constraints  are 

(3.21) 


2      2  > 

ak60  +  Ck6l  =  ° 


ck+i6o2  +  ak6i2  ■  ° 
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where  k  =  0, 1, . . . , (N-l)/2  for  the  a's  and  k  =  0.1,..., N/2  for 
the  c's;  N  is  the  order  of  the  original  driving  point  impedance. 

Performing  the  expansions  above  for  higher-order  augmenting  poly- 
nomials provides  the  results  displayed  in  Table  II. 

The  equations  of  Table  II  can  be  put  into  matrix  format.   Equation 
3.21  becomes 


0    c 


Ck+1   ° 


=  0 


=  0 


(3.22) 


In  general  it  can  be  shown  that 


6fc  Qk  6  >  0 


6C  Rk  6  >  0 


(3.22a) 


where  superscript  t  indicates  a  transpose,  k  =  0,1,..., N/2  and 
j  =  T+l  for  the  first-order  example  (T  is  the  order  of  the  augmenting 
polynomial) . 


Vi 


o 


c.     0 
k 

0     a 


k-1 


"ak   ° 


Ck-1  ° 


'k-1 


0     a 


k-1 


■\    ° 

0     -c 


k-1 


Vl   ° 


0     a 


k-2 


ak-2  ° 


0 


'k-1 


k-2 


k-2 


(3.23) 


24 


k+1 


0  c,  0 

k 


a,  0 

k 

0  c. 


ck    ° 


Vi  ° 


•ck    ° 


k-1 


0        "Vi    ° 

ak-l  °      -ck-l 

0  c 


k-1  k-1     0 


■Ck-1  °  \-2 


=     RV  (3.23) 

-J 


k  k 

Note  that  £.  and  R.  are  dimensioned  (T+l)  x  (T+l) . 

To  illustrate  the  process  so  far  developed  consider  the  dp  impedance 
used  for  the  example  in  chapter  V. 


Z(s)  = 


s  +  l/2s  +  1 


s  +  s  +  1 
which  requires  a  third-order  augmenting  polynomial. 
T  =  3  N  =  2 


ao  =  l 

ax  =  1/2 

CC2  =  1 

e0  =  1 

j 
k 
k 

PjL-1 

=  T  +   1  -  4 
=  0   for  a's 
=0,1   for  c's 

h  =  l 

ao 

= 

axe 

1 

■  a2P0   '  a0P2   = 

-3/2 

ai 

= 

0 

co 

m 

aQ30 

1 

cl 

= 

c^ 

1 
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*" 


-3/2 

0 

-1 

0 

0 

1 

0 

3/2 

-1 

0 

0 

0 

0 

3/2 

0 

0 

24 


0 

0 

-1 

0 

0 

1 

0 

0 

-1 

0 

-3/2 

0 

0 

0 

0 

1 

£- 


1 

0 

3/2 

0 

0 

-3/2 

0 

-1 

3/2 

0 

1 

0 

0 

-1 

0 

0 

4 


1 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

1 

0 

0 

-1 

0 

-3/2 

It  remains  to  solve  for  the  6  coefficients  according  to  eq.  (3.22) 
where 


6  = 


3  11  12 
B.   LINEARIZATION  AND  SOLUTION  OF  CONSTRAINTS  '   ' 

In  this  section  the  non-linear  inequalities  developed  in  the  pre- 
ceding section  are  linearized  using  a  set  of  congruent  transformations 
and  solved.   Inspecting  matrices  0^  and  R,  the  pattern  of  their  formation 

becomes  apparent.   Note  that  they  are  symmetric  for  all  orders.   Evalu- 

k      k 
ating  Q.    and  R.  for  fixed  j=T+l  over  the  indices  k=0, 1 , . . . ,N/2  produces 

a  set  of  constraint  inequalities  as  follow 
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or 


«'  2j  ft  2  0 

6C  R^  6  2  o 
_    _j  _ 


4'  4  i  5  o 


(3.24) 


(3.25) 


6. 


6C  Sf 2   6  2  0 

Inequalities   3.25  represent   N  conditions    involving  T+l  unknowns    in 
(1) 
Where 


■j-1 


As  a  first  step  toward  the  simultaneous  solution  of  the  set  of 


inequalities  above,  consider  the  transforms  for  eq.  3.24 

6  =  M^  x 
6  =  1^  x 


(3.26a) 


k   k  k   k 

where  M.  (N.)  is  the  normalized  modal  matrix  of  Q.  (R.),  i.e.,  the 

k   k 
matrix  of  normalized  eigenvectors  of  £.  (R.).   These  orthogonal 


Since  one  of  the  6.  can  be  chosen  arbitrarily,  there  are  actually  T 
unknowns  in  6 . 
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t  k       k 
transforms  convert  the  quadratic  forms  _6  £•&.(&  R-  $.)  into  a  linear 

combination  of  the  squares  of  the  coordinates  of  the  vector  x  with  no 

cross  products.   Alternately  the  transform  can  be  viewed  as  reducing 

k   k 
the  matrices  (}.  (R.)  to  diagonal  form; 


i.e. 


under  the  transform 


becomes 


6*  £k  6  (3.27) 


6  =  M1^  x  (3.27a) 


x   (M.)  Q.  M.  x  =  x   (M.)    0.  M.  x  =  1  x.   +  1  x.   +  ...  +  1.  ,x.  , 
-   ^-j  *j  -j  -   -   ^-j    * J  - j  -    0  0     11  J-l  j-1 

(3.28) 

k  t     k  -1 
since  a  property  of  an  orthogonal  matrix  is  (M.)   =  (M.)   .   The 

k       k 
l's  in  eq.  3.28  are  the  eigenvalues  of  £.,  and  M.  is  the  orthogonal 

matrix  of  normalized  eigenvectors  which  are  always  available. 

Using  the  transforms  above  in  eq.  3.25  yields 

x   (M.)   Q.  M.  x  =  0 

xfc  (Nk)fc  Rk  Nk  x  =  0 


(3.29) 


Substituting 


Dk  =  (Mk)t  Qk  Mk 

q-J  ^J  *J  "J 

Dk  =  (N*)  Rk  Nk 

r-J  "J  "J  "J 


where  the  D's  are  diagonal  matrices,  3.29  becomes 

x   D.   x  =  0 

-  q-j  - 

x   D.   x  =  0 

-  r-J  - 
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(3.29a) 


(3.30) 


Writing  the  quadratic  forms  as  the  product  of  a  row  and  column  vector 
eq.  3.30  becomes 


k  > 

D.  (row)    z   =  0 

q-j 

D.  (row)    z   =  0 


(3.31) 


where 


z   = 


(3.31a) 


Vi 


k  k  k         k 

and     D.    (row)    (  D.    (row))    indicate      D.    (  D.)   written  as   a   row.      There 
q-j  V-j  q-j      r-j' 

are   N  constraints    on  T+l  variables    in  £. 

The   solution   to   this    linearized   set    of   inequalities   can  be   quite 

2 

elegantly  found  using  Danzig's  simplex  method  modified  to  permit 

optimization  according  to  the  criteria  that  as  many  of  the  individual 
inequalities  as  possible  set  to  zero.   The  simplex  method  initially 
determines  if  the  set  of  inequalities  is  consistent.   For  the  purpose 
of  this  thesis  it  is  recommended  that  one  start  with  T  ■  1.   If  the 
inequalities  are  inconsistent  simply  increase  the  order  of  the  aug- 
menting polynomial  by  one  and  proceed  until  a  feasible  set  is  generated. 

1 
According  to  Guellemin  ,  if  the  even  part  of  Z(s)  has  no  joa-axis  zeros  a 

finite  feasible  set  can  always  be  found. 

A  rigorous  proof  of  this  assertion  is  not  available  in  the  literature 

Next  the  inverse  transformation  of  eq.  3.26a  must  be  accomplished  in  such 

a  way  that  the  resultant  6's  are  equal  (i.e.,  6=6n,  6, =6,,  etc.). 
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It  should  be  noted  that  since  M.  and  NV  are  different,  different 

-J     "J 

6's  can  be  expected  for  the  same  x's.   In  Eq.  3.27  postmultiplying  the 

k   Jc       k    k 
M.  (NV)  by  D   (  D  ),  which  are  arbitrary  diagonal  matrices,  doesn't 
—j  v— j     m—   n— 

alter  the  congruency  of  the  transform  (i.e.,  the  quadratic  form  in  _6  is 

still  transformed  to  a  linear  combination  of  the  squares  of  the  co- 

k   k 

ordinates  of  the  vector  x) .   Since  M.  (N.)  is  the  normalized  modal 

-J  "J 

matrix,  postmultiplying  by  an  arbitrary  diagonal  matrix  is  equivalent 
to  scaling  the  modal  matrix. 

Therefore,  inverse  transforming  eq.  3.26a  in  such  a  way  that  the  6's 

k       k 

are  equal  is  equivalent  to  requiring  that  diagonal  matrices   W  and  W 

be  found  such  that 

k   k      k   k 
M   W  x  =  N   W  x  (3.32) 

-j  m-  -   -j  n-  - 

where  eq.  3.32  is  developed  by  setting  the  two  equations  3.26a, 

postmultiplied  by  the  required  diagonal  matrices,  equal  to  each  other. 

k    k        Jc    k 

For  the  purposes  of  this  thesis  choose  W  =  1   I  and   w  =  p  I , 

n—       —     n-    r  — 

k      k 
where  1  and  p  are  scalars  and  I   is  the  identity  matrix.   With  this 

substitution,  eq.  3.32  becomes 

lk  Mk  x  =  pk  Nk  x  (3.33) 

-J  -      -J  ~ 


or 


^x  =  qk'kNk 
where 


tfT  x  =  q*'   NV  x  (3.33a) 

-J  -  -J  - 


k,k   .k  k 

q    =  1   p 

k  t 
Rewriting  3.33a  by  premultiplying  both  sides  by  (N.)   yields: 

(nV  Mk  x  =  qk,k  x  (3.33b) 
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X 


Expanding  eq.    3.33b   over   the    index  k,(j   =  T+l)   yields 

,*At  w0  0,0  /M0.t  .  1  0,1  /M0.t      N/2  0,N/2 

(N.)      M.    x   =   q    '      x  (N.)      M.x=q'      x  (N.)      M.      x   =   q  x 

-J        "J  _  -  "J        "J   -  ~  "J        "J      ~ 

(N.)   M.  x  =  q  '   x     (N.)   M.  x  =  q  '   x     ...     (N.)   M.   x  =  q  '    x 

yj    - j  -  -      -j    -j  -  -  -j    -j  - 

!  .  .      (3.34) 

,Jl/2.tJ0  N/2,0    ,J/2.t  1     N/2,1         .  N/2.tMN/2     N/2, N/2 

(N\   )M.x  =  q    '  x   (N.   )M.x  =  q     X        (N)M.   x  =  q 

Equations  3.34  are  the  familiar  eigenvalue  problem  in  which  solution 

k  k  k  t   k 

constants  q  '   (the  eigenvalues)  can  always  be  found  for  (N.)  M.  non- 
singular.   Since  equations  3.34  can  always  be  solved  for  the  eigenvalues, 
eq.  3.26a  can  always  be  inverse  transformed  such  that  the  6's  are  equal; 
and,  consequently,  a  vector  _6  can  always  be  found  which  satisfies  each 
of  the  constraints  in  eq.  3.25  simultaneously  once  a  consistent  set  has 
been  generated  and  solved  for  the  vector  x. 

Fortunately  it  is  not  necessary  to  solve  for  6.  explicitly.   It  is 
sufficient  to  realize  that  evaluating  eq.  3.30  at  the  determined  solution 
vector  x  yields  numerical  values: 

x   D.  x  =  k. 

-  q-j  -    Tc 

(3.35) 

t    k     ,   , 
x   D.  x  =  k. 

-  r-j  -    Tc 

t  k 
If  a  particular  scalar  k,  is  zero,  then  the  corresponding  _6  (£.6_ 

t  k 
(6  R.6.)  is  also  zero;  and  if  a  particular  k,  is  positive,  (the  only 

t  k    t  k 
other  possible  condition)  then  the  corresponding  6_  Q._S  (6  R.6.)  are  also 

positive  although  of  different  magnitude.   This  can  be  demonstrated  as 

follows : 

Let 

k    k 
6  =  M .  W  x 

-  -Jr  - 
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which  is 


6  =  lk  Mk  x  (3.35a) 

-J  - 

for 

w"  =  lk  I 

m-  — 

Substituting  eq.  3.35a  in  3.24  yields: 
x'  (M*)'  (lV  fik  Mk  x  =  l0dk)2x02  +  lx  ilh\2  +    ...  +  ij.idVxj.!2 
=  (lk)2  [l0x02  +  l^2  +  ...  +  Ij.^^2]  (3.35b) 

Comparing  eq.  3.35b  with  3.28,  it  is  evident  that  eq.  3.35b  is  eq.  3.28 
multiplied  by  a  positive  constant  which  establishes  the  assertion  above. 

Recalling  that  the  quadratic  forms  are  actually  the  numerator  co- 
efficients of  the  real  part  of  the  augmented  dp  impedance,  by  solving 
eq.  3.31  and  evaluating  3.30  the  corresponding  inequalities  3.24  are 
revealed  as  either  zero  or  positive.   It  will  now  be  demonstrated  that 
this  information  is  sufficient  to  determine  the  optimal  augmenting  and 
augmented  polynomial  to  insure  that  the  real  part  of  Z  (s)  is  positive 
term  by  term. 

C.   OPTIMAL  AUGMENTING  AND  AUGMENTED  POLYNOMIALS 

In  this  section  knowledge  of  the  form  of  the  augmented  driving-point 
impedance  (i.e.,  which  coefficients  are  positive  and  which  zero)  will 
lead  to  the  development  of  the  optimal  augmenting  polynomial  and  con- 
sequently the  optimal  augmented  driving-point  impedance  in  a  straight- 
forward way. 

Consider  an  impedance  function 

m1  +  n1 

Z<S>  =  Z        +  n  (3-36) 

iru  -r  n_ 
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and  an  augmenting  polynomial 


P(s)  =  m  '  +  a  > 


Then 


PCs)        ml     T  "l 
is    the   augmented   polynomial. 


ReZ    (jco)   = 


where 


m. 


"i'V    -   "l""2" 
(m2")2  -    (n2")2 


mm       +  n1n 


S=JC0 


■   =  num.  '   +  n.n,  ' 


ii  . 
"2     ~  "'2"'l 


l2"l 


n,"   =  m,n, '   +  n,m 


ii   . 

li   "  -n 


iu,i 


Substituting  eq.    3.39    in  3.38  yields 


ReZ    (jco)   = 

a 


(■>-■  +  n, ')    (m. 


nl_')    (min^ 


nln2} 


(m^   +  n^)    (m^    -   n^)    (m^ 


n/> 


(3.36a) 


(3.37) 


(3.38) 


(3.39) 


S=JCO 


(3.40) 


ReZ    (jco)   =  Re(P(jco))   ReZ  (jco) 
P(jco) 


(3.41) 


Num(ReZa(ja)))        ^.^   Num(ReZ  (1co)) 


Den(ReZ    (jco)) 


Den(ReZ    (jco)) 


(3.42) 


where   Num(ReZ    (jco))   and  Den(ReZ    (jco))    represent   numerator   and 

3.  3. 

denominator  terms  respectively.   Therefore, 

Num(ReZ  (jco))  =  ReP(jco)  Num(ReZ(jco))  (3.43) 

n 

From  eq.    3.24  we   know  which   coefficients    of  Num(ReZ    (jco))    are   zero 
and  which  are    positive,   and  we  also  know  Num(ReZ(jco))  ,    the    original 
unaugmented  dp   impedance.      A  simple-minded   but   none  the    less   effective 
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way  to  find  ReP(jco)  and  the  unknown  coefficients  of  Num(ReZ  (jco))  is 

8 

to  proceed  by  long  division.   The  scheme  is  best  illustrated  by  a 
numerical  example. 

Given  Num(ReZ(jco))  =  co   -  2co  +  3co  -  4co  +5,  developing  the  form 
of  Num(ReZ  (jco))  using  the  Danzig  method  outlined  above  to  solve  eq. 
3.31  for  z   and  substituting  into  eq.  3.30,  the  following  information 
is  obtained  concerning  Num(ReZ  (jco)) 

Num(ReZ    (jco))   =  co     +0co     +0co +0co  +0co +c  go  +c,, 

where  cr   >  0  and  c,   >  0. 
d  o 

Dividing  Num(ReZ(jco))  into  Num(ReZ  (jco)) 

3. 


+2     +1 


1  -2  +3  -4  +5 


1 1 

1 

+0 
-2 

+0 
+3 

+0 
-4 

+0 
+5 

+C5 

+C5 
+10 

+C6 

+2 
+2 

-3 
-4 

+4 
+6 

-5 
-8 

+1 
+1 

-2 
-2 

+3 
+3 

+  (c  -10) 

-4  D 

+56 

(c  -6) 

(c6-5) 

and  setting  c,.  =  6  and  c,  =  5  the  division  terminates.   The  numerator 
of  the  ReP(jcD)  =  cjd  +  2cd  +1  and  the  numerator  of  ReZ  (jco)  =  00  +6gd  +5. 

3 

A  feature  of  this  method  that  can  be  useful  is  that  the  augmenting  poly- 
nomial ReP(ja>)  always  has  non-negative  coefficients.   It  may  even  be 
Hurwitz .   In  the  Miyata  method  this  added  constraint  is  not  needed;  but 
in  transfer  function  synthesis,  where  augmentation  may  be  called  for, 
the  augmenting  polynomial  is  additionally  required  to  be  Hurwitz. 

It  turns  out  that  if  Z(s)  is  of  second  order  or  less,  it  is  unnecessary 

4 
to  have  exact  knowledge  of  the  form  Num(ReZ  (jco))  since  Poincare  and  Lewis 

have  demonstrated  that  for  Z(s)  quadratic  an  augmented  polynomial  having 

but  three  non-zero  positive  terms  can  invariably  be  generated.   For  Z(s) 

of  greater  than  second  order,  in  general,  exact  knowledge  of  the  form  of 
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Num(ReZ    (jao))    is    required   and   the    preliminary   steps  which  yield    this 

ft 

form  cannot  be  bypassed. 

The  problem  is  now  essentially  solved,  and  a  brief  outline  of  the 
computational  steps  is  in  order. 

First  generate  a  consistent  set  of  non-linear  inequalities  (3.24) 
and  then  employ  the  transform  3.26a  to  set  up  and  solve  the  resultant 
linear  set  of  inequalities  3.31  by  the  Danzig  simplex  method  which 
insures  that  as  many  as  possible  of  the  inequalities  are  zero,  thereby 
minimizing  the  number  of  elements  in  the  subsequent  network  synthesis. 
Armed  with  the  solution  vector  x,  determine  the  form  of  the  solved 
inequalities  3.35  that  yields  which  coefficients  of  the  numerator  of 
the  real  part  of  the  augmented  dp  impedance  are  zero  and  which  are 
positive.   Divide  the  given  Num(ReZ(ja)))  into  Num(ReZ  (jco))  and  solve 

cl 

for  the  unknown  positive  coefficients.   The  process  is  now  complete  and 
the  numerator  of  the  real  part  is  known  (evaluating  the  denominator  is 
straightforward),  and  one  can  proceed  to  the  synthesis  part  of  the 
problem  with  a  polynomial  that  is  positive  term  by  term. 

At  this  point  it  seems  wise  to  demonstrate  the  solution  algorithm 
and  digital  computer  program  developed  to  augment  a  moderately  involved 
driving-point  impedance. 
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IV.   DIGITAL  COMPUTER  PROGRAM 

From  consideration  of  the  computational  scheme  outlined,  it  is 
obvious  that  the  technique  is  admirably  suited  to  implementation  on  a 
digital  computer.   This  is  espacially  true  since  solving  a  set  of  linear 
inequalities  under  the  constraint  that  the  solution  variables  be  greater 
than  zero  is  very  similar  to  the  basic  linear  programming  problem  which 
is  the  subject  of  the  Mathematical  Programming  System   (MPS/360). 

Therefore,  a  highly  developed  and  efficient  system  is  readily 
available  to  solve  the  linearized  inequalities,  and  the  programming 
problem  becomes  one  of  simply  generating  the  inequalities  from  the 
original  driving-point  impedance  and  inputing  the  set  into  the  MPS 
system. 

The  digital  program  required  for  systematic  solution  is  straight- 
forward and,  aside  from  typical  programming  difficulties,  simple.   The 
entire  solution  process  has  not  been  consecutively  programmed  due  to 
interface  restrictions  between  OS/360  and  MPS/360.   The  program  follows 
in  block  diagram  and  Fortran  format. 
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SUBROUTINE  GROAN 
n=l 


<- 


SUBROUTINE  GRUNT 


SUBROUTINE  GROWL 


SUBROUTINE  PUNCH 


n=n+l 


V 
MPS/360 


-> 


Calculate  the  element  values  for  the 
matrix  of  nonlinear  inequalities. 


Form  the  matrix  of  inequalities 


Convert  the  matrix  of  nonlinear 
inequalities  to  a  matrix  of  linear 
inequalities . 


Prepare  an  input  data  deck  for  the  MPS/360 
program  which  solves  the  linear  inequalities 
and  consequently  the  form  of  the  original 
nonlinear  inequalities. 


Output  the  required  order  of  augmentation 
and  form  of  the  numerator  of  the  real  part 
of  the  driving-point  impedance. 
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//  EXEC  FORTCLGP, REGION. G0=100K,TIMF.G0=5 
//FORToSYSIN  DD  * 

DIMENSION  ALF A( 10 ) ,BETA( 10) ,A( 10), C( 10) ,0(10,10,10) 

DIMENSION  R(10,10,10) 

DIMENSION  QR( 10,10)  ,  XX ( 10, 36 ) , YY ( 1 0 ,36) 

1  FORMAT (1 10) 

2  F0RMAT(8F10o2) 
READ  (5,1)  N 
READ  (5,1)  N 
M3=  (M+2)/2 
M4=2*M3 
M5=4*M3 
MM=M+1 
N1=N+1 
N2=N1+M 

READ  (5,2)  (ALFAU  )  ,1  =  1, MM) 

READ  (5,2)  (BETAd  ),I  =  1,MM) 

CALL  GROAN  (M, N2 , M3 , M5 , ALF A, BETA, A, C ) 

D03  N9=1,N1 

CALL  GRUNT  ( M3 ,N9 , A ,C , Q, R ) 

CALL  GROWL  (Q, R, M3 , M4 , N9 ,QR, XX ,YY ) 

CALL  PUNCHY(N2,N9,M4,QR) 

3  CONTINUE 
STOP 
END 

SUBROUTINE  GROAN  ( M ,N2 , M3, M5 , ALFA, BET A, A, C  ) 

DIMENSION  ALFA(IO) , BET A (10) ,A( 10) ,C(10) 

MM=M+2 

DO  1  I=MM,M5 

ALFAU  )=0o0 

1  BETA(I)=0o0 
D02  M1=1,M3 
A(M1  )  =  0o0 
C(M1  )=OoO 
M2=2*M1 

DO  3  N1=1,M2 

I=2*N1 

J=I-1 

K=4*M1-I 

L  =  K+1 

11=1-1 

JJ=J-1 

KK=K-1 

LL=L-1 

IF(K)  4,4,5 

4  K  =  MM 

5  IF(KK)  6,6,7 

6  KK=MM 

7  IF  (LL)  8,8,9 

8  LL=MM 

9  IF  (JJ)  10,10,11 

10  JJ=MM 

11  A(M1)=ALFA( I)*BETA(K)-ALFA( J)*BETA(L)+A(M1» 
C(M1  )=  ALFAU  I  )*BETA(KK)-ALFA(  J  J  )*BETA(  LL  ) +C  (  Ml ) 

3  CONTINUE 

2  CONTINUE 
M6=M3+1 
D012  I=M6,N2 
A( I)=OoO 

12  C(I)=CoO 
RETURN 
END 

SUBROUTINE  GRUNT  ( M3, N9 , A,C ,Q, R ) 

DIMENSION  A ( 10 ),C( 10) ,Q( 10,10, 10), R( 10,10,10) 

DO  1  1=1, M3 

DO  2  J  =  1,NC> 

DO  3  K=1,N9 

Q(I , J,K)=OoC 
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11 

12 

6 
5 
4 


R ( I , J,K) 
CONT  INUE 
CONTINUE 
DO    4    1=1 
L  =  I 
N  =  2*I 
DO    5    J=l 
DO    6    K=J 
Q(I, J,K) 
Q(I,J+1, 
Q(  I, J,K  + 
IF     (J-l) 
Q(I, J  +  l, 
GO    TO    12 
Q( It J+l, 
IFUoEOo 
IF(L-l) 
L  =  L-1 
CONTINUE 
CONTINUE 
DO    15    1= 
L=I 
N  =  2*I 
DO    16 
DO    17 
RUt  J 
R(It 
R(I, 


17 
16 
15 


10 
9 
8 


30 

40 
20 
25 


=  0oC 


,M3 


fN,2 

,N9,4 
=A(L) 
K+1)=C(L) 
2)=-C(L) 

7,7,11 
K+3)=-A(L) 

K+3)=-A(L-l) 
1)     GO    TO    4 
4,5,6 


It  M3 


J  = 
K= 
K) 
J+l, 
J  +  2, 


R(I, J,K+ 
R(I, J+l, 
IFUoEQo 
IF  (L-l) 
L  =  L-1 
CONTINUE 
CONTINUE 
DO  8  1=1 
DO  9  J=l 
M  =  J+1 
DO  10  K= 
Q(I,K, J) 
R(I,K, J) 
CONTINUE 
CONTINUE 
DO  20  1= 
DO  30  J= 
WRITE(6, 
DO  40  J= 
WRITE(6, 
CONTINUE 
FORMAT (/ 
RETURN 
END 


ltN,2 
J,N9,4 
=C(L+1  ) 
K+1)=A(L) 
K+2)=C(L) 
2)=-A(L) 
K+3)=-C(L) 
1)     GO    TO    15 
15,16,17 


,M3 
,N9 

M,N9 

=  Q(  I,J,K) 

=R(  I,J,K) 


1,M3 

1,N9 

25) (0(1, J, K) ,K=1,N9) 

1,N9 

25) (R( I, J,K),K=1,N9) 

/,6F15<,5) 


10 
9 


XX, YY) 
,10) 
, YY(lOf^O) 


SUBRCUTINE    GROWL (Q , R, M3, M4,N9, QR , 

DIMENSION    0(10,10,10) ,R( 10,10 

DIMENSION    QR(lOflO) ,XX(10,40) 

DIMENSION    QC(10) ,RR(10) ,X(40),Y(40) 

DO    12    1=1, M3 

11  =  1 

DO    9    K=1,N9 

DO    10    J=1,K 

QQ(II)=0(I , J,K) 

RR(I I)=R(I , J,K) 

11=11+1 

CONTINUE 

CALL    EIGEX(C0,X,N9,0) 

CALL    EIGEX(RR,Y,N9,0) 

L=l 

I  11  =  2*1-1 

1111=111+1 

DO    15    J=1,N9 


39 


QRU  II  ,  J)=QQ(L) 
QRUIII,  J)=RR(L) 

15  L=J+L+1 
N10=N9**2 

DO    16    J=1,N10 
XX( I,  J)  =  X( J) 

16  YY( I  ,J)=Y( J) 
12    CONTINUE 

DO    40    1=1, MA 
40    WRITE(6,30) (QR(I,J) ,J=1,N9) 

DO    45    1=1, M3 

WRITE (6, 35) (XX (I, J) ,J=1,N10) 
45    WRITE(6,35)  (YY( I, J)  ,J  =  1,N10) 
30    F0RMAT(//,6F15o5) 
35    FORMAT    (//,6F15o5) 

RETURN 

END 


SUBROUTINE  PUNCHY( N2, N9, M4 ,QR ) 
DIMENSION  QR(10, 10) ,QR1( 10,10) 

40  FORMATC  NAME'  ,10X,  'TESTPRQB'  ) 

41  FORMAT( "ROWS' ) 

42  FORMAT  (2X, »N  COST'  ) 

43  F0RMAT(2X,'G  ROW', ID 

44  FORMAT ( 'COLUMNS' ) 

45  FORMAT (4X, 'C0L1' ,6X,'C0ST' ,6X, 'OoO' ) 

46  FORMAT (4X, 'C0L1' ,6X, 'ROW' , I 1 ,6X , F12o 5 ) 

47  FORMAT( 'RHS'  ) 

48  FORMAT  (4X ,' RHS'  ,7X , 'COST' ,6X,  «0o0  •  ) 

49  FORMAT (4X, 'RHS',7X, 'ROW' , 1 1 ,6X , F12o 5 ) 

50  FORMAT( ' ENDATA' ) 

60  FORMAT (4X, 'COL' ,I1,6X, 'COST*  ,6X,  'OoO'  ) 

61  FORMAT (4X, 'COL', II, 6X,  'ROW' , 1 1 , 6X, F12o 5 ) 

62  FORMAT (4X, 'RHS' ,7X,  'ROW', Il,6X,F12o5) 
IF(N9oEQol)    GO    TO    70 

IF    (N9oGTo2  )    GO    TO    30 

WRITE(7,40) 

WRITF(7,41) 

WRITE(7,42) 

WRITE(7,43) (I ,I=1,M4) 

WRITE(7,44) 

WRITE(7,45) 

WRITE (7,46) (I,QR(I ,1) ,I=1,M4) 

WRITE(7,47) 

WRITE(7,48) 

DO    1    1=1, M4 

1  QRK  I,2)=-QR(I,2) 

WRITE  (  7,49  )  (ItORK  1,2)  ,I=1,M4) 
WRITE(7,50) 
GO    TO    70 
30    WRITE(7,40) 
WRITE(7,41 ) 
WRITE(7,^2) 
WRITE(7,43) (I ,I=1,M4) 
WRITE(7,44) 
N8=N9-1 

WRITE(7,60)( 1=1, N8) 
D02    J=1,N8 

2  WRITE  (7,61)  (J,  I,  QRU,  J)  ,I=1,M4) 
WRITE(7,47) 

WRITE(7,48) 
D03    1=1, M4 

3  QRK  I,N9)  =  -QR(I,N9) 

WRITE  (7,  62  )  (I,  QRK  I ,  N9  )  ,  1  =  1 ,  M4  ) 
WRITE(7,50) 
70    RETURN 
END 
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SUBROUTINE     E I GEX ( A, R, N, MV ) 

DIMENSION    A(l)  ,R(1 > 
5    RANGE=loOE-6 

IF(MV-l)  10,25,10 
10    IQ=-N 

DO    20    J=1,N 

IG=IQ+N 

DO    20    1=1, N 

IJ=IQ+I 

R(IJ)=0oC 

IF(I-J)     20,15,20 
15    R( IJ)=loO 
20    CONTINUE 
25    AN0RM=0o0 

DO    35    1=1, N 

DO    35    J=I,N 

IF(I-J)     30,35,30 
30    IA=I+( J*J-J)/2 

ANQRM=ANORM+A(IA)*A( I  A) 
35    CONTINUE 

IF(ANORM)    165,165,40 
40    AN0RM=!o414*SQRT( ANORM) 

ANRMX=ANORM*RANGE/FLOAT(N) 

IND=0 

THR=ANORM 
45  THR=THR/FLO£T(N) 
50  L  =  l 
55  M  =  L+1 
60  MQ=(K*M-M)/2 

LQ=(L*L-L)/2 

LV=L+MQ 
62  IF(  ABS(ACLM) )-THRl  130,65,65 
65  IND=1 

LL=L+LQ 

n^=m+mq 

X=0o5*(A(LL)-A(MM) ) 
68  Y=-A(LM)/  SQRT( A(LM)*A(LM)+X*X) 

IF(X)  70,75,75 
70  Y=-Y 
75  SINX=Y/  SQRT(2oO*(loO+(  SQRT < lo 0-Y*Y ) ) M 

SINX2=SINX*SINX 
78  COSX=  SQRT( Io0-SINX2) 

C0SX2-C0SX*C0SX 

SINCS  =SINX*COSX 

ILQ=N*< l-l) 

IMQ=N*(M-1 ) 

DO  125  1=1, N 

I0=(  1*1-1) /2 

IFU-L)  80,115,80 
80  IF(I-M)  85,115,90 
85  IV=I+MQ 

GO  TO  95 
90  IM=M+IQ 

95  IF(I-L)  100,105,105 
IOC  IL=I+LQ 

GO  TO  110 
105  IL=L+IQ 
110  X  =  A(  IL)*COSX-A(IM)*SINX 

A( IM)=A( IL)*SINX+A(IM)*COSX 

A( IL  )  =  X 
115  IF(MV-l)  12C, 125,120 
120  ILR=ILQ+I 

IMR=IMQ+I 

X  =  R(  ILR)*CCSX-R(IMR)*SINX 

R( IMR)=R(I LR)*SINX+R( IMR)*COSX 

R(ILR)=X 
125  CONTINUE 

X=2oO*A(LM)*SINCS 
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Y=A( LL)*COSX2+A(MM)*SINX2-X 

X  =  A(  LL)*SINX2+A(M^)*COSX2+X 

A(LM)=(A(LL)-A(MM) ) *SI NCS+ A( LM )* (C0SX2-S I NX2 ) 

A(LL)=Y 

A(MM)=X 
130  IF(M-N)  135,140,135 
135  M=M+1 

GG  TO  60 
140  IF(L-(N-1)>  145,150,145 
145  L  =  L+1 

GO  TO  55 
150  IF(IND-l)  160,155,160 
155  IND=0 

GO  TO  50 
160  IF(THR-ANRMX)  165,165,45 
165  RETURN 

END 
//GCoSYSIN  DD  * 
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//JOBLIB    DD   DSNAME=SYSloMPS360LP,DISP=(SHR, PASS) 
//  EXEC  LINPRGG 

//MPSloSYSIN  DD  * 

PROGRAM 

INITIALZ 

MOVE(XPBNAME, 'PBFILE' ) 

MOVE( XDATAt 'TESTPROB' ) 

MOVE(XOBJ, 'COST'  ) 

MOVE(XRHS,  'RHS'  ) 

MVADR(XDONFSt INFEAS) 

MVADR(XDCFEAS»FEAS) 

GCTO( AARRRG) 
INFEAS    MOVE( XPBNAME,  'PBFILE'  ) 

MOVE (XD AT A, 'TESTPROB' ) 

CONVERT( 'SUMMARY') 

GOTO( AAARRG) 
AARRRG         CCNVERT( 'SUMMARY' ) 
AAARRG         BCDOUT 

SETUPCMAX'  ) 

PRIMAL 
FEAS      SOLUTION 

EXIT 

PEND 
//MPS2oSYSIN  DD  * 
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V.      SAMPLE   PROBLEM  AND   SOLUTION 


«(»)  -  %  +  1/2s  *  1 

s     +  s  +    1 


From  computer  program  results 


in         fi         ft  4  9 

NRez    (ico)    =  go       +Ood     +0oo     +   Oco     +c.co     +cc 
a   J  4  5 

4  2 

NRez(jco)      =  co     -   3/2co     +   1 

4        2 
DReZ(jco)      =00     -oo+l 


Dividing  NReZ    (jo))   by  NReZ(joo)   gives: 

3. 


+3/2  +5/4  +3/8 


1  -3/2  +  1 

1 

0 

0 

0 

+c4 

+C5 

1 

-3/2 

1 

3/2 

-1 

3/2 

-9/4 
+5/4 
+5/4 

+3/2 

-3/2 

-15/8 

+3/8 

+3/8 

+5/4 
-5/4 
-9/16 
-11/16 

+3/8 
-3/8 

Choosing  c,  and  c   equal  to  11/16  and  3/8  respectively,  the  division 
terminates  and  the  real-part  numerator  of  the  augmenting  polynomial  is 
seen  to  be 

co6+3/2tD4+5/4oo2+3/8. 

The  real-part  numerator  of  the  augmented  polynomial  has  three  non- 
zero terms  and  reads 

If)       ? 
cd  +11/1603  +3/8. 

The  denominator  is 

DReZ(joo)  x  ReP(joD) 

which   equals 

oo10+1/2co8+3/4cd6+5/8cd4+7/8co2+3/8. 
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Therefore 


10  2 

D  _    ,.   .  a>     +ll/16a>  +3/8 

ReZ    (jod)   = 


a   '  a)10+l/2oo8+3/4a)6+5/8oo4+7/8o)2+3/8 


which   is    positive    for  all  a>,   as    required 
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TABLE    I6 


COMPARISON  OF  THE   NUMBER  OF   CIRCUIT  ELEMENTS    REQUIRED  FOR  SINGLE-N 
AND    1/2 -N  SPLIT  MIYATA  AND  MODIFIED  BOTT-DUFFIN  METHODS 


Single-n 

Half-n 

Modified 

Split 

n 

r 

lc 

Split 

n 
r 

Bott-Duffin 

n 

!k 

n 

r 

1 

i 

2 

1 

2 

1 

2 

2 

4 

3 

4 

3 
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3 

3 

9 

4 

6 

4 

7 

5 

4 

16 

5 

10 

5 

15 

7 

5 

25 

6 

14 

6 

19 

11 

6 

36 

7 

17 

7 

35 

15 

8 

64 

9 

25 

9 

75 

31 

10 

100 

11 

35 

11 

155 

63 

12 

144 

13 

44 

13 

315 

127 

14 

196 

15 

52 

15 

635 

255 

16 

256 

17 

62 

17 

1275 

511 

18 

324 

19 

74 

19 

2555 

1023 

20 

400 

21 

86 

21 

5115 

2047 

30 

900 

31 

139 

31 

163,835 

65,535 

n1      =  number   of   reactive   components 
n        =  number   of  resistive   components 
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